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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 

RESEARCH MEMORANDUM 

THE INTERRELATED EFFECTS OF ENGINE ROTOR MOMENTUM AND FLIG [T 
REGIME ON THE DYNAMIC STABILITY OF HIGH-PERFORMANCE AIRCRAFsi 

By Warren J. North 



SUMMARY 

A five-degree-of -freedom analysis of aircraft maneuvering stability 
showed that gyroscopic moments due to engine rotor momentum produced 
considerable aircraft rolling asymmetry for certain values of stability 
derivatives and flight conditions. 

Because of its powerful effect on inertia-coupling and effective dihe- 
dral, the most important consideration was the initial angle of attack. 
Critical reductions in dynamic stability were caused by variations of 
stability derivatives with angle of attack. 

At high altitudes and Mach numbers the effective dihedral due to 
swept wings and large dorsal fin caused oscillations which were either 
divergent or lightly damp ed. The periods of oscillations encountered at 
high altitude are of the same order of magnitude as pilot reaction time. 

It appears, therefore, that artificial damping must be incorporated in 
many high-performance aircraft. 


INTRODUCTION 

The stability of rolling aircraft has become a critical design and 
operational consideration with the advent of supersonic, high-altitude 
aircraft. As predicted in reference 1, concentration of weight in the 
long fuselage and reduced aerodynamic restoring moments at high altitude 
have caused inertia- coupled rolling divergence in some cases or have in- 
creased aircraft yaw and pitch during roll to a point where structural 
redesign has been required. Flight experiences with roll-coupling are 
discussed in references 2 and 3. Analog computer studies of transonic 
rolling stability are given in references 3 and 4. 

•^The information presented herein was offered as a thesis in partial fulfill- 
ment of the requirements for the degree of Master of Aeronautical Engineering, 
Princeton University, Princeton, New Jersey, December, 1956. 

IFIED 
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Since the trend in supersonic aircraft design is toward thin wings 
of short span, most of the aircraft and fuel weight is concentrated 
along the fuselage. Therefore, during a rapid rolling maneuver the air- 
plane tends to rotate about its principal longitudinal axis rather than 
the wind axis. The first 90° of roll translates the initial angle of 
attack into a sideslip, the magnitude of sideslip depending on Mach num- 
ber, directional stability, altitude, and rate of roll. The cross- 
coupling aerodynamic, inertia, and elastic terms then complicate the 
maneuver and require a theoretical analog which involves all the air- 
plane degrees of freedom. 

The ratio of engine thrust to airplane weight is increasing as air- 
plane design speeds increase; turbojet engines designed for very high 
altitudes will need large cross-sectional areas to handle sufficient 
mass flow; vertical- takeoff aircraft need large thrust-to-weight ratios 
and will operate at airspeeds corresponding to near-zero aerodynamic 
forces and moments. In light of these trends, the engine rotor gyro- 
scopic moment will become a more important dynamic stability parameter, 
since, in general, the rotor momentum will increase as engines increase 
in size and thrust. 

In references 1, 5, and 6, the assumption of steady rolling was 
made in order to linearize the equations of motion. Reference 5 extends 
the theory of reference 1 to show the effect of engine rotor momentum on 
steady-rolling stability. Reference 6 also extends the theory of refer- 
ence 1 and demonstrates how the transients in angle of attack and side- 
slip can be approximated using the steady-rolling assumption. During a 
supersonic tactical maneuver the rolling velocity will probably be a 
continual transient since there will be a few instances when it will be 
necessary to roll rapidly to angles greater than 180°. This report 
analyzes the transient airplane response during a half-roll maneuver 
for a wide range of Mach number and shows the effects of variations in 
engine rotor momentum, Mach number, altitude, static stability, and 
load factor. 

The basic airplane chosen for this theoretical investigation was a 
supersonic interceptor with highly swept delta wing and tail surfaces. 


SYMBOIS 


CD 



wing span, ft 

drag coefficient, drag/q Q S 

lift coefficient, lift/q c S 

rolling-moment coefficient, rolling moment/qQSb 
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C m 

pitching -moment coefficient, M/q Q £>c 

c n 

yawing-moment coefficient, N/q^Sb 

C Y 

lateral-force coefficient, lateral force/q^ 

c 

wing mean aerodynamic chord, ft 

4394 
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V* 

forces with respect to Y and Z axes, lb 

s 

acceleration due to gravity, ft/sec 2 


lx > ~^-Y> I Z moments of inertia about X, Y, and Z principal axes ; slug- 



sq ft 

Ix e 

M 

§ , 

•° H 

• — 1 

moment of inertia of engine rotating parts about X axis, 
slug-sq ft 

stabilizer deflection, deg 

o l w 

incidence of symmetrical-airfoil wing, deg 

L 

rolling moment, ft- lb 

M 

pitching moment, ft-lb 

m 

mass of airplane, W/g, slug 

N 

yawing moment, ft-lb 

P 

rolling velocity, radian/sec 

p c 

critical roll rate, radian/sec 

Pe 

engine rotor angular velocity, radian/sec 

<1 

pitching velocity, radian/sec 

<Jo 

dynamic pressure, lb/sq ft 

r 

yawing velocity, radian/sec 

m 

s 

wing area, sq ft 

u,v,w 

components of velocity V along X, Y, Z principal axes, 
ft/sec 
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V velocity, ft/sec 

W airplane weight, lb 

X,Y,Z airplane principal axes 
a angle of attack of principal X axis, deg 

a + e, aerodynamic angle of attack, deg 
0 angle of sideslip, deg 

5 a ,S r aileron and rudder deflection, respectively, deg 

e inclination of body axis above positive X axis, deg 

0 Euler elevation angle of positive X axis, deg 

4> Euler roll angle, deg 

i f Euler yaw angle, deg 

Superscript: 

derivative with respect to time 
Stability derivatives: 
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PROCEDURE 

The motion of a rigid airplane can be completely described by 
eight differential equations which represent the six degrees of freedom 
and involve conversion from Euler angular velocities to angular veloci- 
ties of principal axes. The initial time history of a rolling maneuver 
can be represented satisfactorily by seven equations if the airplane 
velocity is assumed nearly constant during these first few seconds of 
the maneuver. The Mach number will decrease slightly due to the com- 
pound effect of greater induced drag and decreased engine thrust associ- 
ated with inlet flow distortion. However, this Mach number change can 
be neglected if the airplane is not operating at or just above the tran- 
sonic range where there is a rapid variation of aerodynamic coefficients 
with Mach number. The remaining seven equations are derived in appendix 
A. The equations are written in terms of the principal axes; therefore, 
the product of inertia terms do not appear in the moment equations. 

| When the equations of motion are referred to principal axes, air- 

plane drag might be expected to be an important stability consideration, 
especially at the low lift-drag ratios associated with supersonic flight. 

I Calculations were made with and without the drag terms in the equations 

for lift and side force. The only significant effect of drag was a 
slight damping of the extremely large angle-of -attack excursions. 

| 

The static stability derivatives for a typical supersonic inter- 
ceptor were obtained from wind tunnel tests or Mach number extrapolations 
thereof. Damping derivative components due to the tail surfaces were 
calculated from static derivatives, whereas wing and fuselage damping 
components were obtained from theory. Aeroelastic motions were not 
treated as additional degrees of freedom but in this analysis were con- 
sidered as fixed corrections to the stability derivatives. Since wind 
tunnel balances obtain some forces and moments in terms of wind or sta- 
bility axes, the aerodynamic derivatives were transferred to the princi- 
pal axes where necessary. Four of the derivatives were considered func- 
tions of angle of attack. The stability derivatives and drag coefficients 
for the basic airplane are shown in table I. Airplane dimensional and 
mass constants are shown in table II. 

The amount of engine rotor momentum used in this analysis was con- 
sidered to be the maximum for the assumed aircraft size and performance 
capabilities. 


Machine Computation 

In the interest of accuracy and in order to include many nonlinear 
derivatives and functions, the IBM 650 digital computer was used. A 
picture of the dual machine setup at the NACA Lewis laboratory is shown 
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in figure 1. The Runge-Kutta method of numerical integration was se- 
lected. With an integration step size of 0.04 second, 2 hours of machine 
time were required to compute and card-punch the time history of a 6- 
second maneuver. The cards were tabulated and also read into an auto- 
matic curve-plotting machine. 


Control Input 

The aileron deflection was programmed independently as shown in 
appendix B. Deflection was dependent on roll angle to the extent that 
roll position determined the point at which aileron neutralization began. 
For each reference run (zero engine rotor momentum) a trial-and- error 
procedure was used whereby duration of aileron deflection and amount of 
reverse aileron were chosen in order to roll to approximately 180° as 
rapidly as possible. 

During the initial stages of the investigation an attempt was made 
to maintain roll angle with corrective aileron after the desired roll 
was attained. This corrective aileron deflection was attempted 
using various combinations of roll velocity, roll acceleration, and roll 
position as controlling parameters. A pilot reaction time of 0.24 second 
was assumed. The results showed that in most cases this pilot analog was 
sufficiently out of phase with the maneuver so as to aggravate the error 
in roll angle. In addition, under extreme angle of attack, effective 
aileron reversal existed, and the situation was again aggravated when 
corrective aileron was programmed. Consequently, in order to provide a 
comparison during left and right rolls in which engine rotor momentum 
was the only asymmetry, the ailerons were programmed neutral subsequent 
to the reversal deflection. 

Several rolling maneuvers were continued to angles greater t.han 180° 
in an attempt to compare divergence tendencies with the steady-roll 
theory of reference 1. 

At the higher values of dynamic pressure, aileron deflection may be 
limited by large aileron forces and hinge moments. However, for the 
purpose of this report, momentary full deflection was programmed at all 
flight conditions. 


Although the initial value of stabilizer position was dependent on 
flight condition, the stabilizer position was maintained constant 
throughout each maneuver. Zero rudder deflection was maintained for all 
maneuvers . 
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RESULTS AND DISCUSSION 

Airplane dynamic behavior during rolling maneuvers is complicated 
by variations of flight parameters and interrelations thereof. There- 
fore, when discussing the effect of a single parameter such as engine 
rotor momentum on aircraft dynamic stability, the corresponding flight 
conditions and stability parameters must be specified. 

This report attempts to show first the effects of rotor momentum 
associated with changes of altitude, Mach number, and initial load 
factor. Then, for cases where the effect of rotor momentum was most 
significant, some of the stability parameters were varied in order to 
determine specific combinations of flight conditions and stability param- 
eters for which the effect of rotor momentum was most critical. 


Aeroelastic effects are concealed in Mach number and altitude ef- 
fects. Figure 2 shows that equal dynamic pressures of 1300 pounds per 
square foot exist at Mach numbers of 1.7 and 3.5 and altitudes of 30,000 
and 60,000 feet, respectively. At this dynamic pressure, reductions of 
10 to 60 percent in static stability and control effectiveness are caused 
by aeroelastic deformation. The effect of aeroelasticity on the sta- 
bility derivatives can be seen in table I. 


Altitude Effects 


Half-roll maneuvers were calculated for altitudes of 30,000 and 
60,000 feet at a common Mach number of 1.5. The initial load factor 
was 1 g. The maneuver time histories at 30,000 feet are shown for left 
roll without rotor momentum and for left and right rolls with rotor mo- 
mentum in figures 3, 4, and 5, respectively. There was little varia- 
tion in maximum angle of attack or sideslip for variations in rotor mo- 

1 ° 

mentum. The maximum excursion in angle of attack was 3- , and maximum 

c* 


sideslip was -2~ . At 60,000 feet (fig® 6) a half-roll maneuver caused 

a maximum excursion in angle of attack of 8.8° and sideslip of -6.8°. 
Left roll (fig. 7) with engine momentum resulted in slightly higher val- 
ues of angle of attack and sideslip, while right roll (fig. 8) was ac- 
companied by slightly lower values. Similar engine effect trends would 
be expected from the steady- roll theory of reference 5. Although the 
engine momentum had a small effect on the maneuver at both altitudes, 
the variation was somewhat greater at the higher altitude. Since the 
magnitude of maximum roll rate was 4.5 radians per second at the lower 
altitude as compared with 1.6 at the higher altitude, it appears that 
the roll rate had lesser effect on angle of attack and sideslip than 
did altitude. The altitude effect was largely due to reduction of dy- 
namic pressure which caused an increase in initial inclination of the 
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principal axis from 3° to 9.3°. The reduction of roll rate with in- 
creased altitude is due to reduction in dynamic pressure, reduction in 
aileron effectiveness at high angles of attack, and strong dihedral ef- 
fect (positive rolling moment due to negative sideslip). In fact, at 
the higher altitude (fig. 6) the magnitude of rolling velocity de- 
creased to zero at 2 seconds although full aileron was applied. Note 
also that angle of attack had increased to 13° and angle of sideslip was 
negative. Since dihedral effect also increases with angle of attack for 
a swept midwing configuration, increasing angle of attack has a compound 
effect on reduction of rolling velocity. 


Mach Number Effects 

The effect of supersonic Mach number on a half-roll maneuver can 
be determined by comparing figures 6 (Mach number, 1.5) and 9 (Mach 
number, 3.5). Both figures represent 60,000 feet and level flight. At 
the higher Mach number the aircraft responded much more rapidly to ai- 
leron input since angle of attack was small and dynamic pressure was 
high. The small excursions in angles of attack and sideslip at the 
higher Mach number are due to small initial inclination of the principal 
axis and large restoring moments. Although the excursions of angles of 
attack and sideslip were smaller at Mach number 3.5, the increments of 
normal and lateral accelerations were greater because of a five-fold in- 
crease in dynamic pressure. At 1 second, during the high-speed maneuver, 
the angle of attack and corresponding airplane normal acceleration became 
negative. The negative acceleration sensed by the pilot would be aug- 
mented by a negative increment because of his position above the roll axis. 
From the pilot's viewpoint this negative acceleration, coupled with the 
higher normal and lateral accelerations, would make the high-speed 
maneuver less tolerable. 

It is interesting to note that the period of oscillation in side- 
slip i- s essentially the same at both Mach numbers. This is probably due 
to the reduction in directional stability at the higher Mach number. 

The inclusion of engine rotor momentum had no significant effect on the 
maneuver at Mach number 3.5. 


Initial Load Factor 

From the standpoint of inertia-coupling, variation of initial load 
factor is synonymous with variation of the principal axis inclination. 
Figures 10, 9, and 11 show calculated half-roll maneuvers at a Mach 
number of 3.5 and an altitude of 60,000 feet for initial load factors 
of 0, 1, and 2, respectively. The corresponding initial inclinations 
of the principal axis were 1°, 4°, and 7.3°. At zero g initial condi- 
tion (fig. 10), the maximum roll velocity exceeded 5 radians per second 
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magnitude, but primarily because the principal axis was essentially 
alined with the flight path during the maneuver, there was little vari- 
ation in angles of attack and sideslip as the airplane rolled. The ab- 
sence of yawing moment due to aileron deflection aided in permitting 
this extremely stable rolling condition. At 1 g initial condition 
(fig. 9) it has been observed previously that there was a moderate os- 
cillation, but the variation in angles of attack and sideslip were 
within allowable limits. At 2 g initial condition (fig. 11) large os- 
cillations are noted in sideslip and roll, which were triggered by 
inertia-coupling and maintained by strong effective dihedral coupling. 
Angle of attack increased beyond the allowable value corresponding to 
7 g maximum airplane load factor. Figures 12 and 13 present left- and 
and right-roll maneuvers at 2 g's with the engine rotor momentum in- 
cluded. It can be noted that variations in angles of attack and side- 
slip are increased slightly in left rolls and retarded in right rolls. 

Increased load factor was not as critical at lower dynamic pres- 
sures. At low pressures an increase in load factor required large angles 
of attack. At high angles of attack the effective dihedral counteracted 
much of the aileron rolling moment so that high roll rates became more 
difficult to attain. 


Effective Dihedral 

It has been previously noted that the basic airplane configuration 
exhibited strong positive dihedral effect. Figure 14 shows a half-roll 
maneuver with zero effective dihedral at a Mach number of 1.5 and an 
altitude of 60,000 feet. Zero engine rotor momentum was assumed during 
this maneuver. Zero effective dihedral might be realized with a low 
straight-wing installation incorporating a ventral fin in addition to 
the conventional dorsal fin. Negative wing dihedral would also tend to 
offset the yaw- induced rolling moment due to dorsal fin or angle of 
attack. As noted in figure 14, the initial angle of attack was trans- 
lated into negative sideslip during the first 90° of roll. Positive 
directional stability would normally decrease sideslip, but the com- 
bined effect of inertia- coupling and decreasing directional stability 
caused the motion to diverge in sideslip, reaching -30° at 3.9 seconds. 
The directional stability decreased because of increasing angle of at- 
tack. Figures 15 and 16 show left and right rolls with the inclusion of 
rotor momentum. If the airplane equations of motion are analyzed, it can 
be noted that the following engine rotor momentum terms are added to the 
pitch and yaw equations, respectively: 

^gPe I X e Pe 

- r ’ ~TT~ q 


CONFIDENTIAL 


4394 


CU-2 back 


NACA RM E57B05 


CONFIDENTIAL 


11 


At time subsequent to 2.4 seconds (fig. 15) the positive pitching ve- 
locity causes an engine-induced positive yawing moment which accelerated 
the divergence in sideslip. The sideslip reached a magnitude of 30° in 
3.5 seconds in this left-roll maneuver. During the right roll (fig. 16) 
different phase relations are found between pitching and yawing veloc- 
ities . As a consequence , the rotor couple differs from that for the 
left roll and tends to stabilize the maneuver. Although not as pro- 
nounced, similar dynamic tendencies were found at a Mach number of 0.85 
and an altitude of 30,000 feet when the dihedral effect was zero. In 
comparing figure 14 with figure 6, the omission of effective dihedral 
uncouples the aerodynamic forces such that rolling velocity is not im- 
peded and also has a destabilizing effect, since sideslip is not bled 
off in roll. Aircraft stabilization due to positive effective dihedral 
is shown in figure 11. This run represents a Mach number of 3.5 for 
which directional stability was very low and decreased to zero at an 
angle of attack of 1 ® , Although the angle of attack during a majority 
of the run was greater than 7°, the airplane did not diverge in side- 
slip but developed a Dutch-roll-type oscillation with coupling in roll 
and sideslip. 


Static Stability C, 


m > 
a 


CL ^ • and C„ 


An important drag consideration for aircraft which operate over a 
wide range of Mach numbers is the large trim drag due to high longitudi- 
nal stability at low supersonic Mach numbers. In-flight reduction of 
longitudinal stability would therefore be desirable. One way of reduc- 
ing longitudinal stability in flight would be by fuel transfer. Figure 
17 shows the airplane response to aileron roll in level flight at a Mach 
number of 1.5 and an altitude of 60,000 feet with longitudinal stability 
reduced by a factor of 10. By comparing figures 17 and 6 the initial 
effect of the reduced stability is seen to reduce the angle of attack 
to lower values during the initial roll transient. The lower angle of 
attack permitted greater aileron effectiveness, which created a faster 
roll and a sharper recovery. The subsequent coupling in roll and yaw 
caused a severe slowly damped oscillation. Angle of attack slowly di- 
verged during the maneuver reaching 33° at 6 seconds. Note that the 
period of the small superimposed oscillation in angle of attack was one- 
half that in yaw and roll. The periods of the oscillations decreased 
considerably during this brief maneuver. The slow divergence in angle 
of attack was, in part, due to the increase in pitching moment with 
sideslip C m ^ 2 . Figure 18 shows that the omission of C m ^ 2 stabilizes 

the maneuver. 


Reference 1 indicates that for steady-rolling maneuvers the critical 
roll rates above which the aircraft would exhibit divergence in angles of 
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attack or sideslip can be calculated. These critical roll rates are 
the natural frequencies of the independent longitudinal and lateral 
modes. However, reference 1 shows that an aircraft should never diverge 
in a steady-rolling maneuver if the critical roll rates in pitch and 
yaw are identical. The permissible difference between the two critical 
rates is shown to increase with damping in pitch and yaw. The critical 
roll rates for the basic configuration are shown in table III. Since di- 
rectional stability was dependent on angle of attack, the value of direc- 
tional stability used in these critical roll calculations corresponded to 
the steady-state angle of attack prior to aileron deflection. For the 
original stability levels used in figures 7 and 8, level flight, a Mach 
number of 1.5, and an altitude of 60,000 feet, the corresponding critical 
roll rates in pitch are -1.68 and 1.80 radians per second, respectively, 
for left and right rolls with rotor momentum. The critical roll rates in 
pitch and yaw are nearly equal. Although peak roll rates slightly exceed 
the calculated critical values for the reference airplane, average roll 
rates were considerably lower and the convergent oscillation might have 
been anticipated from steady-rolling divergence theory. The critical roll 
rates in pitch were reduced to -0.53 and 0.57 radian per second when lon- 
gitudinal stability was reduced by a factor of 10 (figs. 19 and 20). The 
maximum roll rates during the initial roll exceed the critical rates by 
a factor of four, yet this portion of the maneuver did not exhibit a 
divergence. The slow divergence in angle of attack occurred afted the 
aileron was neutralized. 

Figures 21, 22, and 23 show airplane behavior when both directional 
and longitudinal stability were reduced to low values. The pitching and 
yawing critical roll velocities were chosen equal. Again the roll rate 
during the initial portion of the run exceeded the critical rates by a 
factor of four. Recovery from the roll caused divergent oscillations in 
yaw and pitch and a divergent increase in angle of attack. The diver- 
gence was more severe than that associated with unequal static stabili- 
ties. The effects of rotor momentum were small but similar to previously 
mentioned trends. 


Roll Velocity 

Figures 24, 25, and 26 show the time histories of rolling maneuvers 
in which aileron deflection was maintained until high rates of roll were 
realized. These runs were made at a Mach number of 3.5 and an altitude 
of 60,000 feet. The initial values of normal acceleration were the same 
as those of the previous section Initial Load Factor. It is noted that 
for zero g (figs. 10 and 24) there is little variation in angles of 
attack and sideslip for roll rates of 5 and 10 radians per second. At 
1 g initial acceleration (figs. 9 and 25) the variation in angles of 
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attack and sideslip increased slightly as maximum roll rate increased 
from 5.2 to 9.5 radians per second. At 2 g initial acceleration (figs. 
11 and 26) the amplitude of angles of attack and sideslip increased to 
divergence as roll velocity was increased from a magnitude of 4.3 to 
8.4 radians per second but then became stable as roll velocity was in- 
creased to 12.5. Apparently the extremely high rolling velocity had a 
spin- stabilizing effect on the airplane. It is also noted in figure 26 
that the largest variations in angles of attack and sideslip occurred 
after the aileron was neutralized. In every case mentioned previously 
the maximum roll rate exceeded the critical roll rate. Although the 
calculated critical values may be used as an approximate criteria for 
divergence, it appears that, for the transient- type rolling maneuver 
considered herein, initial inclination of the principal longitudinal 
axis will be the important consideration. 


Comparison with Simplified Longitudinal- and Lateral-Mode Theory 

independent longitudinal and lateral modes were calculated by 
the linear small-disturbance theory. All stability derivatives were 
assumed constant and were based on angle of attack at the beginning of 
the maneuver. The independent modes can be compared with the computed 
motions subsequent to the half roll when the ailerons were returned to 
neutral . 

The longitudinal equations consisted of only the lift and pitch 
equations since airplane velocity was assumed constant. The resulting 
quadratic equation might be expected to give two imaginary roots repre- 
senting a short-period oscillation. However, for all the flight condi- 
tions itemized in table I, the longitudinal quadratic gave two real 
roots representing two aperiodic modes , one with a fast convergence and 
one with a weak convergence. The five- degree- of -freedom analysis for 
Mach number 3.5 (fig. ll) showed a slow oscillation in angle of attack 
with a faster oscillation superimposed. The oscillations were probably 
due to the inertia- and aerodynamic -coupling caused by rolling velocity. 
These disturbing forces were neglected in the linear solution for the 
longitudinal mode. 

Calculations of the lateral mode were made for two Mach number con- 
ditions at an altitude of 60,000 feet. The linearized side force, roll, 
and yaw equations were combined to form a quartic from which four roots 
were extracted. For each of the two cases investigated the quartic 
yielded two real roots and two imaginary roots. In each case, the real 
roots represented a slow divergence and a slow convergence. At Mach 
number 1.5 the imaginary roots indicate an oscillation with a 4- second 

period which damps to half-amplitude in 2-| seconds. At Mach number 3.5 
the lateral-mode theory predicts a period of 8 seconds which is damped 
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to half -amplitude in 1 second. At Mach numbers of 1.5 and 3.5 (figs. 6 
and 11) the five-degree-of -freedom analysis shows periods which are 
smaller by factors of 3 and 8, respectively. The actual damping of the 
oscillations varies during the maneuver considered. The divergences 
in some of the motions would not be predicted by linear theory, because 
the destabilizing parameters such as initial inclination of principal 
axis, nonlinear derivatives, and inertia-coupling are not considered 
in the linearized theory. 

Because of the many simplifications which are made in the linear 
theory, it is not surprising that there is poor agreement between this 
and the machine-computed results. This disagreement stresses the need 
for at least f ive-degree-of-freedom calculations when computing the 
transient dynamics of high-altitude, high-performance aircraft. 


CONCLUDING REMARKS 

A theoretical analysis of dynamic stability for high-performance 
aircraft showed that several combinations of flight parameters caused 
divergence or large oscillations in angles of attack and sideslip as a 
result of a half-roll maneuver. The initial inclination of the princi- 
pal longitudinal axis was the most critical consideration. There were 
no roll-induced oscillations when the initial load factor was zero, that 
is, when the principal longitudinal axis was nearly alined with the 
flight path. 

In one supersonic right-roll maneuver at low dynamic pressure the 
addition of engine rotor momentum caused gyroscopic moments sufficient 
to stabilize a divergence. However, inclusion of the rotor momentum 
during left roll aggravated the divergence. The corresponding airplane 
configuration possessed zero effective dihedral. This is an important 
stability consideration because low values of effective dihedral may be 
prevalent in future supersonic airplanes. The addition of a ventral 
fin, negative wing dihedral, or the incorporation of a straight wing 
would reduce airplane effective dihedral. The use of a ventral fin may 
be required to offset the decrease in supersonic directional stability 
associated with the dorsal fin at increasing angles of attack. 

For most configurations and flight regimes investigated the inclu- 
sion of engine rotor momentum caused larger oscillations during left 
roll than during right roll. This asymmetric roll response will re- 
quire different pilot techniques dependent on roll direction. 

An important consideration when computing dynamic stability is the 
variation of stability derivatives with angle of attack. Increasing 
angle of attack can cause reductions or reversals in aileron effective- 
ness, yaw due to roll, and directional stability. It also has a 
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powerful effect on swept- wing effective dihedral. These angle-of- 
attack effects can cause dynamic instability and variations in frequency 
of the oscillatory motions. 

When using electronic computers to compute aircraft maneuvering 
stability, divergences and oscillations should be assessed in terms of 
pilot controllability. As in the cases of the classical phugoid and 
spiral modes, which were often divergent but controllable, some of the 
supersonic roll-coupled instabilities may be controllable. For the air- 
plane configuration considered herein, the fast lateral oscillations, 
primarily of the Dutch-roll type, would be difficult to control. How- 
ever, the large angles of attack usually resulted from slower pitching 
motions which would be easier to control. Assuming that pitch control 
was feasible, the angle of attack should not become excessive and the 
amplitudes of the angle-of-attack-dependent lateral motions should be 
reduced. Since the critical motions occurred after recovery from the 
initial roll, pilot judgment and control would not be affected by the 
radial acceleration associated with rapid rolling, which occurs during 
the initial portion of the maneuver. 

From the standpoint of pilot technique, maneuvering divergences and 
oscillations can be decreased by performing the roll portion of the 
maneuver at reduced load factor. This technique may retard slightly the 
desired change in flight path direction. However, in the case of an in- 
terceptor, the absence of subsequent roll-coupled oscillations would 
permit a more stable gun or missile-launching platform. The above pro- 
cedure would not be applicable if the principal longitudinal axis were 
considerably below the body axis. 

The incorporation of large positive effective dihedral caused 
severe roll-and-yaw coupled oscillations which were either divergent or 
lightly damped. Since the periods of the oscillations were of the same 
order of magnitude as pilot reaction time, it appears that artificial 
damping must be incorporated on many high-altitude, high-speed aircraft. 

Poor correlation of results with simplified linear longitudinal - 
and lateral -mode theory stresses the need for at least five-degree -of - 
freedom calculations when computing the transient dynamics of high- 
altitude, high-performance aircraft. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, February 11, 1957 
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APPENDIX A 


DERIVATION OF EQUATIONS 
Principal Axis System 

The relation of the airplane axes with the velocity and gravity 
vectors is shown in the following sketch: 


X 




(a) 


CONFIDENTIAL 


4394 


NACA RM E57B05 


CONFIDENTIAL 


17 


Equations of Motion 

Velocity components: 

u = V cos a cos p 
v = V sin p 
w = V sin a 

Lift equation: 


■"inertia 


= -m(w - uq + vp) 


■J- = VP cos p 
w = Vd cos a 


' ^gravity 


aero 


= -mVa cos a + mVq cos a cos P - mVp sin P 
= mg cos 9 cos <{> 


= -C^S cos a " ^ ( lo s cos a " CdIo 13 sin 

"b 


a 


Equating ZF z = 0 and solving for d: 

4l ■ * 00sP ' p + f - ar - ar % - % g gH 

(Al) 

Side force: 

F„ = -m(v - wp + ur) 

1 inertia 

= -mVp cos p + mVp sin a - mVr cos a cos p 

F y = mg sin <t> cos 9 

•^gravity 


^aero = q o^°rB P + ^o S gy C Y r r " C D<lo 3 sin 3 


o _ sin a 
P = P , Q - r cos 


cos P 

V % p 


, g sin <t> cos 9 
a *V cos 3 


c ^o 3 sin P 
D mV cos P 


q Q Sb 


mV cos p gmV 2 Y r cos 3 


(A2) 
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Rolling moment: 
^inertia = 

^aero = 

/. P = 

Pitching moment 
^inertia = 

^aero ~ 

/. q = 

Yawing moment: 
Ninertia = 

N 

aero 
/. r = 


-plx + q.r(ly - Iz) 

q^Sb 2 q^Sb 2 

^o SbC l 6a 5 a + ~ 2V C J p P + ~ 2V C l / + % SbC lp P 

Iy - Iz qoSb cuSb 2 qoSb 2 q„Sb 

j- <l r + i Cl 8 a Sa + 2VI X Cz p P + 2VI X Cl r T + "TT Cl 3 P 


(A3) 


-qI Y + pr(l z - I x ) - I Xe P e r 

q^Sc 2 q^Sc 2 

<3o^ C m a L + it + ~W~ C m <1 + ~zy— C mA + ^ 0 ^^ 

a 1 ^ q a (3^ 

P Z - P X Ix e Pe 9oSc „ „ . 

__ pr _ __ r + — Cma a L + — ^it + 




^Sc . q^Sc 


2VIy Cm q q + 2VIy Cm d a + Iy ^P 213 


(A4) 


-rlZ + (I x - Iy)P ( l + IXgPe 1 ! 

q^Sb 2 q^Sb 2 

a SbC 3 + — C r + — C p + q SbC 6 + qSbC 6 

np H T 2V n r 2V n 5 a H n 5 r 


I x - I Y % Pe q Sb q^Sb 2 

T z pq + 4 + ~1~ Cn P P + 2VI Z Cn r r + 


q^Sb 2 qgSb q^Sb 

W Cn p P + ~W Cn 5a 5a + — Cn br &T 


(A5) 
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Conversion of Angular Velocities from Principal to Euler Axes 

The order of rotation of the axis system is shown in the following 
sketch: 


cr> 

to 


M 

o 

cd 

y 

to 

!=) 

o 



In order to keep the gravity force properly oriented during a maneu- 
ver, the relation between the velocities of the Eulerian and principal 
body axes must be determined. If an airplane is assumed initially in 
level flight and then is displaced arbitrarily, the displacement can be 


considered a succession of three rotational velocities: 

(l) Rotation 

* 

about Z axis 

(2) Rotation 

• 

e 

about Y axis 

(3) Rotation 

i 

about X axis 
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The resultant body angular velocity vectors p, q, and r are then 
considered orthogonal. The body rotational velocities can be written 
in terms of the nonorthogonal Eulerian rotational velocities as follows: 

• • 

p = <t> - \|r sin 9 

• • 

q = 9 cos 4> + \|r sin <J) cos 9 

• • 

r = \|r cos 4> cos G - G sin 4> 

A simultaneous solution of the above three equations yields the Eulerian 
angular velocities in terms of the body angular velocities: 

4> = p + q. tan 9 sin + r cos <t> tan 9 (A6) 

% 

9 = q cos <J> - r sin <}> (A7) 

* , cos 4> , sin 4* 

t - r tt + q - 

cos 9 cos 9 

Only equations (A6) and (A7) are required to describe the change in 
Eulerian angles which affect the gravitational force, gravity being in- 
dependent of airplane directional sense . 

In addition to the five differential equations of motion, equa- 
tions (A6) and (A7) were also included in the set which was solved on 
the computer. 
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APPENDIX B 


AILERON SCHEDULE 


The rate of total aileron deflection was chosen as 100° per second 
with .±28° as maximum deflection: 

& a = ±100, 0 (depending on 4> , & a ) (Bl) 

A sample left-roll aileron schedule is illustrated in the following 
sketch: 



Time 
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TABLE I. - STABILITY DERIVATIVES AND DRAG COEFFICIENTS 


Altitude, ft 

30,000 

60,000 

Mach number 

0.85 

1.5 

1.5 

3.5 


per deg 

0.056 


0.0419 

0.0508 


0.0225 

° L i t 

per deg 

0.0087 


0.00508 

0.0060 


0.00286 

Si 

per deg 

-0.00645 


-0.00846 

-0.0124 


-0.001536 

C m it 

per deg 

-0.01413 


-0.00916 

-0.0109 


-0.00527 

S 

per rad 

-4.17 


-2.78 

-3.22 


-1.45 

Sx 

per rad 

-1.69 


-0.498 

-0.715 


-0.557 

c v 

per deg^ 

0.00033 


0.00033 

0.00033 


0.00025 

% 

per deg 

-0.0202 


-0.01665 

-0.01955 


-0.01176 


per rad 

1.03 


1.08 

1.08 


0.32 

CZ P 

per rad 

-0.208 


-0.110 

-0.1732 


-0.0585 

c l 

L r 

per rad 

0.185 


0.117 

0.117 


0.061 

°n r 

per rad 

-0.885 


-0.925 

-0.925 


-0.350 

Cn 5 a 

per deg 

0 


0 

0 


0 

% 

per deg 

-0.0029 - 
0.00008 


-0.00275 - 
0.00015 

-0.00275 - 
0.00015 


-0.0008 - 
0.00002 

X 

per deg 

0.001042 - 
0.302X10” 5 

a L 

0.000472 - 
0.12X10" 5 

0.00085 - 
0.216X10” 5 

a L 

0.0003 - 
0.294X10” 6 a£ 

S 

per rad 

0.045 - 
0 . 0138 ccj ^ 


0.0303 - 
0.0128 cx^ 

0.0356 - 
0.0151 


0.0292 - 
0.0063 

% 

per deg 

0.00635 - 
0.207X10” 4 
(a L > 0) 

a L 

0.00459 - 
0.141X10” 4 a? 
(a L > 0) 

0.0054 - 
0.166X10” 4 
(a L > 0) 

a L 

0.001 - 
0.00015 

C D 

0.0155 + 
0.000488 a£ 

0.024 + 
0.000415 a l 

0.024 + 
0.000582 

0.020 + 
0.000358 a£ 
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TABLE II. - AIRPLANE CONSTANTS 


S, sq ft 

401 

m, slug 

1,087 

b, ft 

35.8 

c, ft 

15 

e, deg 

-1 

i w > deg 

0 

Iy, slug-sq ft 

11,000 

Iy, slug-sq ft 

345,000 

Ig, slug-sq ft 

344,000 

IXgPg^ slug-sq ft/sec 

40,000 


TABLE III. 


Yaw: p c = ± 


Pitch: P c = ± 


- CRITICAL ROLL RATES 


1 



1 

/■ c m a ( lo S ^\ 2 J X e ^e 
\ I Z " X X/ + 2 ( J Z " 


Mach 

number 

Altitude, 

ftxio -3 

Initial 

load 

factor 

Critical roll rate, 
radian/sec 

Yaw 

Pitch 

Right 

Left 

Right 

Left 

0.85 

30 

1 

2.18 

-2.07 

1.52 

-1.40 

.85 

30 

2 

1.83 

-1.71 

1.52 

-1.40 

1.5 

30 

1 

3.37 

-3.25 

3.00 

-2.88 

1.5 

30 

2 

3.29 

-3.18 

3.00 

-2.88 

1.5 

60 

1 ' 

1.64 

-1.52 

1.80 

-1.68 

1.5 

60 

2 

.65 

-.53 

1.80 

-1.68 

3.5 

60 

1 

1.37 

-1.25 

1.50 

-1.38 

3.5 

60 

2 

.53 

-.41 

1.50 

-1.38 
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Figure 3. - Left-roll maneuver at Mach number of 1.5, altitude of 30,000 feet, and 
level flight. 
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Figure 4. - Left-roll maneuver at Mach number of 1.5, altitude of 30,000 
feet, and level flight. Rotor momentum included. 
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Figure 5. - Right-roll maneuver at Mach number of 1.5, 
feet, and level flight. Rotor momentum included. 
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Figure 6. - Left-roll maneuver at Mach number of 1.5, altitude of 60,000 
feet, and level flight. 
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Figure 7. - Left-roll maneuver at Mach number of 1.5, altitude of 60,000 
feet, and level flight. Rotor momentum included. 
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Figure 8. - Right-roll maneuver at Mach number of 1.5, altitude of 60,000 
feet, and level flight. Rotor momentum included. 
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Time, sec 

Figure 9. - Left-roll maneuver at Mach number of 3.5, altitude of 60,000 
feet, and level flight. 
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Figure 10. - Left-roll maneuver at Mach number of 3.5, altitude of 60,000 
feet, and zero g pushover. 
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Figure 11. - Left-roll maneuver at Mach number of 3.5, altitude of 60,000 
feet, and 2 g pull-up. 
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Figure 12. - Left-roll maneuver at Mach number of 3.5, altitude of 60,000 
feet, and 2 g pull-up. Rotor momentum included. 
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Figure 13. - Right-roll maneuver at Mach number of 3.5, altitude of 60,000 
feet, and 2 g pull-up. Rotor momentum included. 
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Figure 14. - Left-roll maneuver at Mach number of 1.5, altitude of 60,000 

feet, and level flight. C. =0. 
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Figure 15. - Left-roll maneuver at Mach number of 1.5, altitude 
of 60,000 feet, and level flight. C = 0; rotor momentum 
included. 
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Figure 16. - Eight-roll maneuver at Mach number of 1.5, altitude of 60,000 feet, and level 
flight. C 7 = 0; rotor momentum included. 
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Figure 17. - Left-roll maneuver at Mach number of 1.5, altitude of 60,000 feet, 
and level flight. = -0.00124. 
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Figure 18. - Left-roll maneuver at Mach number of 1.5, altitude of 60.000 feet 
and level flight. = -0.00124; C m ^ 2 =0. ' ' 
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Figure 19. - Left-roll maneuver at Mach number of 1.5, altitude 
and level flight. C m ■ -0.00124; rotor momentum included. 
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Figure 20. - Right-roll maneuver at Mach number of 1.5, altitude of 60,000 feet, 
and level flight. C m = -0.00124; rotor momentum included. 
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Figure 21. - Left-roll maneuver at Mach number of 1.5, altitude of 60,000 

feet, and level flight. C m = -0.00124; C n = 0.000536 - 0.166X10 -4 a 

ni a 
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Figure 22. - Left-roll maneuver at Mach number of 1.5, altitude of 60,000 feet, 

and level flight. C m = -0.00124; C n = 0.000536 - 0.166X10 -4 a^; rotor 

ct £3 

momentum included. 
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Figure 23. - Right-roll maneuver at Mach number of 1.5, altitude of 60,000 

feet, and level flight. C = -0.00124; C n = 0.000536 - 0.166X10 4 a?; 

a “6 Jj 

rotor momentum included. 
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Figure 24. - Left-roll maneuver at Mach number of 3.5, altitude of 
60,000 feet, zero g pushover, and large roll rate. 
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Figure 25. - Left-roll maneuver at Mach number of 5.5, altitude of 
60,000 feet, level flight, and large roll rate. 
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(a) Large roll rate. 


Figure 26. - Left-roll maneuver at Mach number of 3.5, altitude of 
60,000 feet, and 2 g pull-up. 
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(b) Extreme roll rate. 

Figure 26. - Concluded. Left-roll maneuver at Mach number of 
3.5, altitude of 60,000 feet, and 2 g pull-up. 
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